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FOREWORD 

The  work  described  in  tlus  e  xport  was  petfonned  in  the  Space  and  Surface  Systems  Divimon  of 
the  Strateipe  and  Space  Systems  Department  at  the  request  of  the  Cruise  Mianle  Weapon  Systems 
Division  (LIO)  of  the  Strike  Systems  Department.  A  description  of  the  analysis  and  software  developed 
to  find  the  intersection  of  two  convex  polygons  is  pven. 

The  intersecti<m  of  convex  polygons  was  required  by  Silnlle  Tallant  of  Lll  during  her  study  to 
determine  operating  areas  for  strikes  against  multiple  targets.  This  algorithm  is  used  in  a  strike 
analysis  tool  that  she  developed  called  STANT.  The  author  has  benefited  frmn  numerous  discuasions 
with  her  on  the  application  of  the  algorithm  and  is  <q>preciative  of  her  thorough  review  of  this  report. 
Mrs.  Tallant  also  produced  the  drawings  that  appear  in  this  document. 
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ABSTRACT 

An  algorithm  is  pven  that  finda  the  inteiaection  of  two  convex  pcdygona.  It  ia  coded  in  Fiwtran 
for  the  IBM  PC  deaktop  computer.  The  program  ia  robuat  and  laat.  It  haa  been  uaed  aucceaafhlly  in 
targeting  applications  that  require  a  r^id  determination  of  the  common  intersection  of  more  than  100 
cmivex  polygona  each  apecified  by  more  than  150  verticea. 


ui/iv 


NSWCDD/TR-93/345 


CONTENTS 

Page 

l.  INTRODUCTION . 1 

n.  SUBROUTINE  XINT:  MASTER  ROUTINE  (PART  1)  . 2 

m.  FUNCTION  ITR:  ALGORITHM  A . 4 

IV.  SUBROUTINE  FST:  ALGORITHM  B  (with  SORT) . 5 

SUBROUTINE  SORT  -  CALLED  BY  FST . 6 

V.  SUBROUTINE  SOLV:  ALGORITHM  C . 6 

VI.  SUBROUTINE  XINT  (PART  2)  . 7 

FUNCTION  INT  -  CALLED  BY  XINT . 8 

REFERENCES . 9 

APPENDIX -EXAMPLES  USING  FINTX,  BASED  ON  INTSEC . A-1 

DISTRIBUTION . (1) 


NSWCDD/TRr#3/345 


L  INTRODUCTION 

This  report  describes  an  algorithm,  INTSEC,  that  determines  the  r^ion  intersecticm,  INTX, 
of  taro  convex  closed  polygons.  We  shall  refer  to  the  mtire  PC  Fmtran  code  for  INTSEC  by  the  name 
FINTX.  The  code  is  robust;  it  will  always  find  INTX,  except  in  cases  where  the  single  precision 
arithmetic  of  the  PC  cannot  resolve  <»  distinguish  between  different  points  (a  double  (wecision  version 
of  the  code  that  provides  greater  resdution  is  also  available).  FINTX  is  also  fisst.  In  a  recent 
application,  using  an  IBM  compatible  486-M  DX2  deskti^  computer,  FINTX  found  the  conunon  inter¬ 
section  of  172  polygons  each  with  180  vertices  in  6.0  seconds.  Examples  using  FINTX  are  given  in  the 
appendix. 

INTSEC  has  important  applications  in  computer  gra4>hk8,  computer  chip  design,  and  targeting 
studies.  For  example  in  targeting,  INTSEC  is  usdul  in  generating  an  operating  area  against  several 
targets.  Some  papers  studied  on  the  intersection  problem  are  given  in  references  [1],  [3],  [4],  [5],  [6]. 
None  of  these  peters  gave  sufficient  detui  of  the  actual  implementation  of  their  algorithms  to  evaluate 
their  speed.  For  example,  an  algorithm  may  be  carried  oat  in  a  small  number  of  steps  but  if  each  step 
is  expensive  time-wise,  as  in  computing  arctangents,  its  efficiency  is  reduced.  No  remarks  were  found 
concerning  robustness  of  their  algorithms. 

The  only  requirement  of  the  two  given  polygons,  in  addition  to  being  closed  and  convex,  is 
that  they  be  positively  oriented  (PO).  A  rimple  polygon  P  is  positively  oriented  (PO)  if  its  vertices  are 
ordered  with  the  interior  of  P  on  the  Jgfl  as  the  boundary  of  P  is  traversed  in  the  direction  of  increasing 
indices  k  of  the  vertices  Pj^,  k  =  1,2,  N.  If  j[j^  is  replaced  by  right  then  P  is  negatively  oriented  (NO). 
The  two  polygons  are  denoted  by  XY  and  UV,  where  XY  is  spedfied  by  its  vertices  Sj  with  coordinates 
(xj,  yj),  i  =  l,2,...,  NX.  Similarly,  UV  is  specified  by  its  vertices  Wj  with  coordinates  (uj,  Vj), 
j  =  1,2,...,  NU.  Since  the  polygons  are  closed  (xj,  y,)  =  (Xffxy  Vnx)  (“i»  =  (“NU’  ^Nv)- 

our  procedure  we  also  include  {Xf,x  +  v  ysx -I- 1)  =  (*»’  72)  (“iVl/  + 1’  ''NU -f-l)  =  (“2»  '^2)-  ^ 

remainder  of  this  report  this  fact  will  not  be  referred  to  explicitly. 

Note  that  INTX  is  also  convex  and  is  completely  determined  by  its  vertices.  It  can  have  at 
most  KK  of  them  with  KK  <  NXl  -i-  NUl,  where 

NX1  =  NX-1>3,  NU1  =  NU-1>3. 

The  vertices  will  be  ordered  so  that  INTX  is  PO  with  the  first  vertex  taken  as  the  <me  that  has  the 
smallest  ordinate.  If  nune  than  one  such  vertex  exists,  then  the  one  of  that  set  that  has  the  minimum 
abscissa  is  taken  as  the  first  vertex. 

INTSEX7  is  made  up  of  three  basic  algorithms  (A,  B,  C)  and  some  auxiliary  algorithms. 
Algorithm  A  establishes  if  a  given  vertex  w  of  UV  is  either  inside  XY,  outside  XY,  or  on  its  boundary, 
dXY.  It  also  determines  if  a  given  vertex  s  of  XY  is  inside  UV,  outside  UV,  or  (m  its  boundary,  dUV. 
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Algorithm  B  finds  the  first  intersection  of  dXY  and  dUV  that  has  not  been  found  by  A.  Algorithm  C 
finds  all  the  remaining  vertices  of  INTX  not  found  by  A.  A  vertex  of  INTX,  which  u  not  a  vertex  of 
either  XY  or  UV,  is  found  by  C  at  the  intersection  of  a  line  segment  (edge)  of  XY  and  a  line  segment 
(edge)  of  UV,  where  a  line  segment  (edge)  is  specified  by  the  coordinates  of  its  end  points. 

For  easy  reference,  we  list  here  the  names  of  the  routines  associated  with  the  algorithms 
mentioned  above.  They  will  be  discussed  in  the  order  listed. 


ROUTINE  ALGORITHM  SECTION 


XINT  (PART  1) 

ITR 

FST 

SORT 

SOLV 

XINT  (PART  2) 
INT 


Master  routine 
ALGORITHM  A 
ALGORITHM  B 
Called  by  FST 
ALGORITHM  C 
Master  routine 
Called  by  XINT 


Described  in  Section  II 
Described  in  Section  III 
Described  in  Section  IV 
Described  in  Section  IV 
Described  in  Section  V 
Described  in  Section  VI 
Described  in  Section  VI 


No  proofs  are  given  in  this  paper. 


n.  SUBROUTINE  XINT.  MASTER  ROUTINE  (PART  1) 


The  first  part  of  the  master  or  executive  routine  XINT  is  described  here.  Before  proceeding, 
some  notation  is  introduced.  A  given  polygon  P  may  be  defined  by  the  sequence  of  its  vertices  {pj^}, 
k  =  1,2, ...,  N,  where  P  is  generated  by  taking  the  pj^  in  increasing  order  of  the  subscripts.  The 
statements  q  is  in  P,  q  is  contained  in  P,  or  q  e  P  means  q  is  located  inside  P  or  on  its  boundary,  dP. 
As  noted  earlier  we  refer  to  XY  by  its  vertices  {s,-},  i  =  1,2,...,  NX  or  simply  by  {s}.  Similarly,  UV  is 
defined  by  {w},  or  {w^},  j=:l,2,...,  NU.  Note  that  XY  has  NXl  =  NX  —  1  distinct  vertices  and  UV  has 
NUl  =  NU  —  1  such  points.  No  preliminary  processing  of  XY  or  UV  is  necessary. 


XINT  begins  by  calling  the  Fortran  function  ITR(u,  v,  NXl,  X,  Y),  which  is  based  on  a  part  of 
algorithm  A,  to  find  if  the  vertex  w  s  (u,  v)  of  UV  is  contained  in  XY,  where  X  and  Y  are  arrays 
containing  the  x  and  y  coordinates,  respectively,  of  {s}.  ITR(u,  v,  NXl,  X,  Y)  can  take  one  of  three 
values  that  is  stored  in  Lw.  If  Lw  =  1,  then  w  is  in  XY  but  it  does  not  coincide  with  a  vertex  of  XY.  If 
Lw  =  0,  then  w  coincides  with  a  vertex  of  XY.  If  Lw  =  - 1,  then  w  is  not  in  XY.  XINT  determines  Lw 
for  each  w  in  UV.  If  ITR  finds  NUl  vertices  of  UV  in  XY,  then  INTX  has  been  found,  namely 
INTX  =  UV.  If  this  is  not  the  case,  XINT  calb  nR(x,  y,  NUl,  U,  V)  to  find  if  the  vertex  s  s  (x,  y)  of 
XY  b  contained  in  UV,  where  U  and  V  contain  the  coordinates  of  the  {w}.  The  result  from 
ITR(x,  y,  NUl,  U,  V)  b  stored  in  Lz.  Ls  takes  one  of  the  three  values  (1,  0  ,  or  - 1)  with  meanings 
analogous  to  those  for  Lw.  XINT  determines  Lz  for  each  z  in  XY.  If  NXl  vertices  of  XY  are  found  in 
UV  by  ITR,  then  INTX  =  XY.  The  abscbsas  and  ordinates  of  the  vertices  belonging  to  INTX  found  by 
ITR  are  stored  in  the  arrays  WW  and  ZZ,  respectively.  In  order  to  avoid  duplications,  if  Lz  =  0  then 
the  coordinates  of  the  vertex  are  not  saved,  since  they  have  already  been  found  by  Lw  =  0,  w  s  z  and 
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■torad  in  WW  and  ZZ.  At  this  point,  we  Mnime  tl«^  WW  and  ZZ  each  contain  K1  elonents,  with 
0  <  K1  <  NXl  4-  NUl.  We  nae  WW(K)  and  ZZ(K)  to  denote  the  Kth  elenMnt  theee  arrays. 

XINT  now  calls  FST,  which  is  based  on  algorithm  B.  FST  searches  for  the  first  intersection 
point  c  of  dXY  and  dUV  that  is  not  a  vertex  of  either  of  the  polygrms.  Let  SiSj  denote  the  directed  line 
segment  frmn  end  point  S|  to  end  pcwt  s^.  The  search  for  c  b^ins  with  directed  edges  S}S2  and  WjWj. 
If  their  intersection  produces  c  an  exit  occurs,  otherwise  the  search  continues  by  kxking  for  a  smallest  j 
such  that  2  <  j  <  NUl,  intersects  Mfy-  V  no  c  is  found,  the  procedure  is  rq>eated  with  j  =  1 

and  with  SiSj  replaced  by  s,-s,-  ^  |,  i  =  2.  The  entire  procedure  is  continued  by  incremoiting  j  through 
its  range  for  each  i  <  NXl  until  c  is  found  for  the  smallest  values  of  j  and  i.  If  no  c  is  found  then  INTX 
has  sero  area  and  an  exit  is  made.  If  c  exists  for  smne  smallest  i  s  n  and  j  »  m,  this  implies  that  there 
is  an  edge  of  XY,  and  an  edge  of  UV,  that  intersect  at  c.  Recall  that  c  is  not  a 

vertex  of  XY  or  UV.  The  vertices  of  XY  and  UV  are  then  reordered  by  the  auxiliary  routine  SORT 
such  that  the  sequences  {>}  and  {w}  are  rotated  as  shown: 

{*n»  *n  +  V  •••«  *1»  ***»  •!»}  "*  {*!»  *3*  **•»  *Nx) 

{’'m.  *m  +  l.  •••• 

More  details  on  FST  are  given  in  Section  IV. 

After  finding  c  XINT  calls  SOLV,  which  m  based  on  algorithm  C.  SOLV  finds  the  remaining 
intersection  points  that  make  up  INTX  by  moving  around  dXY  and  dW  and  systematically  finding 
the  misring  points.  The  detmls  of  this  search  by  SOLV  are  given  in  Section  V.  After  all  vertices  of 
INTX  have  been  found,  the  routine  HULL,  in  [2],  is  used  to  reorder  the  pmnts  as  described  in 
Section  I  and  to  store  their  coordinates  in  new  output  arrays  W  and  Z.  Then  a  fin*l  search  is  made  in 
W  and  Z  for  any  successive  duplicate  points;  if  any  exist,  tmly  one  set  of  such  points  is  retained.  The 
arrays  W  and  Z  h<dd  the  comdinates  of  the  wdered  vertices  of  INTX. 

An  example  of  how  INTSEC  operates  is  shown  by  Figure  1. 


FIGURE  1.  AN  EXAMPLE,  INTX  FOUND  BY  INTSEC 


3 


NSWCDD/TR-93/545 


Points  a,  b  are  found  in  that  order  by  ITR.  Then  point  c  is  found  by  FST.  Note  that  since  c  is  found 
by  FST,  the  points  of  UV  are  reordered  by  SORT  such  that 

Wj  -*  Wj,  Wj  -♦  Wj,  -*  W3,  Wj  -*  W4,  Wj  -♦  Wj, 

or  as  stated  above 

(Wj,  Wj,  W4  ,  Wj,  Wj)  ^  (W|,  W2,  Wj  ,  W4,  Wj). 

Then  SOLV  is  called  to  find  the  remaining  points  d,e,f  INTX.  HULL  is  then  used  to  reorder  the 
points  of  INTX  yielding  INTX  =  {a,  c,  b,  d,  e,  f,  a). 

The  remainder  of  XINT  will  be  described  in  Section  VI.  The  order  with  which  the  array 
elements  of  XY  and  UV  are  presented  to  SOLV  for  finding  the  intersectiem  pt^ts  of  the  edges  of  XY 
and  UV  that  have  not  been  obtained  by  ITR  is  discussed. 


m.  FUNCTION  ITR:  ALGORITHM  A 

Given  a  PO  polygon  P  specified  by  its  vertices  {pj,  pj,  ...,  p^-,  ...,  Pj^}  and  a  point  q,  ITR 
determines  if  q  e  P  with  q  ^  Pj  for  each  j,  if  q  =  p^-  for  some  j,  or  if  q  is  outside  P.  From  ITR  a 
parameter  Lq  is  assigned  a  value  1,  0,  - 1,  accordingly.  The  evaluation  of  Lq  by  algorithm  A  depends 
strongly  on  the  following: 

Let  A  denote  a  triangle.  Then  S(A),  or  S  for  short,  has  the  properties  that  |S |  is  twice  the  area 
of  A  and  can  be  given  in  terms  of  the  Cartesian  coordinates  ((,  17)  of  the  vertices  of  A.  Specifying  A  by 

{Pi.  P3.  P3t  Pi)  P,  =  »Ki)] . 

S  =  [^(2)  -  ((!)]  (q(3)  -  m  +  [€(3)  -  m  [nil)  -  ni2)] .  (1) 

If  A  is  PO  then  S  >  0,  and  if  A  is  NO  then  S  <  0;  S  =  0  imidies  the  vertices  are  colinear. 

Again,  we  have  for  a  PO  polygon  P  and  a  point  q,  with  N1  =  N  - 1: 

a)  Lq  =  1.  q  (  P,  but  not  at  a  vertex  of  P. 

b)  Lq  =  0.  q  =  p^-  for  some  j  =  1,  2, ...,  Nl. 

c)  Lq  =  - 1.  q  /  P. 

Let  q  =  (3c,  y).  Algorithm  A,  using  ITR,  begins  by  checking  to  see  if  q  =  pj.  If  so,  then  Lq  =  0  and  an 
exit  from  ITR  occurs.  Otherwise,  it  keeps  the  line  segment  qp|  fixed  and  proceeds  counterclockwise 
around  P  looking  at  the  sequence  of  triangles  A^-  =  {q,  P|,  p^,  q},  j  «  2, 3, ...,  Nl.  Starting  with  j  s  2,  the 
quantity 

S(A,)  =  («(l)-3f]  [q(j)-yl  +  Ua)-Xl  lF-q(l)l  (2) 

is  evaluated  to  determine  the  orientation  of  A^.  For  rim|dicity  in  notation  let  S(A^)  =  S^-.  If  for  each  j 
S^-  >  0,  then  q  is  not  in  P  and  Lq  =  - 1.  If  there  exists  an  integer  k  such  that  2  <  k  <  Nl  and  S^  <  0, 
then  an  additional  new  triangle.  A,  is  considered.  It  is  defined  by  its  sequence  of  vertices  (q,  pj^  _  j,  p|^, 
q)  and  its  orientation  is  determined  by 

S(A)  =  §*_!  =  [«k-l)-3fl  [i7(k)-Fj  +  U(k)-fJ  (y-.j(k-l)].  (3) 
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Ufiiig  these  concq>te  we  eanuMrise  the  leeulte  of  the  verioue  poeeibilitiee. 
If  q  «  Pj,  then  Lq  «  0. 

If  >  0  for  each  j,  then  Lq  =  - 1. 

If  Sj^  <  0  end  _  i  <  0i  then  Lq  «  - 1. 

If  <  0  and  _  i  >  0|  then  Lq  >  1. 

If  a  0  and  _  |  >  0,  that  Lq  « 1. 

If  Sjt  a  0  and  |  «  0,  then  Lq  «  0. 


IV.  SUBROUTINE  FST:  ALGORITHM  B  (with  SORT) 


Subroutine  FST  finds  the  first  intosection  c  ot  an  edge  XY  with  an  edge  of  UV  that  is  not 
at  an  end  point  either  edge.  The  procedure  bepns  by  looking  for  the  first  intersectUHi  of  the  ith  edge 
*(*i  It  >  of  ot  edges  w^w^  ^  j,  for  increasing  j,  j  s  1,  2, NUl  of  UV.  If  one  is  not  found, 

then  i  is  incronented  by  1  and  the  process  is  repeated.  If  for  i  >  NXl  no  intersecticm  has  been  found, 
then  INTX  has  area  0  and  an  exit  is  made  fimn  XINT.  Thus,  let  c  =  ((,  q)  denote  the  intersection 
point  of  the  ith  edge  of  XY  with  the  jth  edge  of  UV,  where  the  end  points  of  the  ith  edge  have 
coordinates  [x(i),  y(i)]  and  [x(i-f  1),  y(i  +  l)],  and  the  end  points  of  the  jth  edge  have  coordinates 
[u(j),  v(j)]  and  [u(j  + 1),  v(j  -f  1)] .  Then  the  equations  to  be  satisfied  ate 

Dxq~Dy(sB 


where 


Du  q  —  Dv  (sC, 


Dx  =  x(i  -t- 1)  —  x(i) 

Du  =  uQ  +  l)-u(j) 

B  =  y(i)x(i  +  l)-x(i)y(i  +  l) 

^  _  C  Dx  -  B  Du 

TliEL  — 


Dy  =  y(i+i)-y(») 

Dv  =  v(j  +  l)-v(j) 

C  =  v(j)  u(j-)-l)-u())  v(j  +  l) 
CDy-BDv 

’  = - DEE - 


Now,  let 


DEL  =  Dy  Du  —  Dx  Dv. 


(4) 


Ts|DxDu|  +  iDyDv|. 

If  |DEL|  < TE,  E  =  e/4  =  1.25*10~^,  then  the  two  edges  under  omsideratimi  are  numerically  parallel 

and  cannot  yield  c.  Hence,  assume  that  |DEL|  >  TE.  We  check  to  see  if  ((,  q)  is  contained  in  the 

rectangle  R(e)  specified  by  the  inequalities 

Xmn  -  e  |Xmn|  <  (  <  Xmx  ■+■  e  |Xmx| 

Ymn  — c  |Ymnj  < q< Ymx  +  e  |Ymxj , 

where 

Xmn  s  max[min(x(i),  x(i  -f  1)),  min(u(j),  u()  -t- 1))] 

Xmx  =  min[max(x(i),  x(i-t- 1)),  max(u()),  u(j  + 1))] 

Ymn  s  max[min(y(i),  y(i-f  1)),  min(vQ),  v(j  +  1))J 
Ymx  s  min[max(y(i),  y(i  + 1)),  max(v(j),  v(i  + 1))] . 
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The  value  of  e  is  chosen  as  a  unall  oS  the  smallest  positive  single  |»edn(m  number  in  the  IBM 

PC  fm  which  1  +  e  >  1.  Figure  2  shows  lt(0). 


FIGURE  2.  RECTANGLE  R(0) 

SUBROUTINE  SORT-CALLED  BY  XNIT 

If  ((,  I})  is  in  R(e)  then  it  is  accepted  as  the  first  intersection  pmnt,  provided  it  does  not 
coincide  with  an  end  point  of  dther  edge  which  can  occur  in  spite  of  the  iact  that  such  pmnts  have 
already  been  found  by  ITR.  Assuming  c  has  been  found,  the  polygras  are  reordered  by  the  SORT 
routine.  This  routine  is  most  eadly  described  by  sim|dy  listing  its  few  lines  of  Fortran  code,  which  is 
done  below  for  the  XY  polygon.  Suppose  Sj^j^  ^  j  denotes  the  XY  edge  of  the  intersection.  SORT  sU»cs 
x(k)  in  x(l)  as  the  first  element  of  the  sequence,  x(k+l)  in  x(2),  etc.,  as  described  earlier.  Similarly,  the 
y(i)  are  reordered  in  the  same  way.  Hence  if  X  and  Y  are  the  arrays  to  be  reordered,  then  SORT 
requites  as  input  NXl,  X,  Y,  and  k.  The  algorithm  is  given  by 


SUBROUTINE  SORT(NXl,  X,  Y,  k) 

(DIMENSION  STATEMENT  FOR  X,  Y,  XI,  Yl) 

D05  M  =  l,k-1 

X1(M)  =  X(M)  !  XI  AND  Yl  ARE  TEMPORARY  STORAGE  ARRAYS 

5  Y1(M)  =  Y(M)  !  WITH  THE  SAME  DIMENSIONS  AS  X  AND  Y. 

L  =  0 

DO  10  M  =  k,  NXl 
L  =  L+1 
X(L)  =  X(M) 

10  Y(L)=Y(M) 

N2  =  NXl-k  +  l 
DO  15  M  =  l,  k-1 
X(N2  +  M)«X1(M) 

15  Y(N2  +  M)  =  Y1(M) 

X(NXl  +  l)sX(l) 

Y(NX1  +  1)«Y(1) 

X(NXl-t-2)sX(2) 

Y(NX1+2)»Y(2) 

END 
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V.  SUBROUTINE  SOLV:  ALGORITHM  C 

After  using  FST,  XINT  calls  SOLV  in  order  to  find  the  remaining  intersection  points  of  INTX. 
SOLV  takes  as  input  the  coordinates  of  the  end  points  of  a  directed  edge  of  XY  and  the  end  points  a 
directed  edge  of  UV  and  determines  if  the  two  lines  cross.  The  coordinates  of  the  intersection  point  are 
given  as  output  and  a  parameter  MO  is  assigned  one  of  the  output  values:  1,  2,  3.  If  MO  =  1,  then  the 
lines  do  not  cross.  If  MO  =  2,  then  the  lines  cross  inside  their  end  points  and  the  intersection  point  is 
accepted  as  a  new  point  of  INTX.  If  MO  =  3,  then  the  two  edges  overlap  and  the  crossing  point 
coincides  with  at  least  one  of  the  four  end  points  that  has  already  been  found  by  ITR.  The  procedure 
to  determine  the  crossing  point  is  the  same  as  described  for  FST  and  the  same  equations  hold. 

Of  course  once  MO  returns  a  value,  XINT  must  determine  how  to  proceed  and  not  miss  any 
remaining  intersection  points.  If  MO  =  3,  then  we  have  found  a  previous  intersection  point,  so  we  treat 
it  as  a  crossing  point,  as  if  MO  =  2,  but  do  not  store  it  in  our  INTX  arrays  WW  and  ZZ.  Hence,  it  is 
sufficient  to  have  a  way  of  proceeding  in  the  two  different  situations  MO  =  1  and  MO  =  2.  The 
procedures  are  described  in  Section  VI  with  the  use  of  the  function  INT. 

VL  SUBROUTINE  XINT  (  PART  2) 

In  this  section  we  describe  how  the  remaining  points  of  INTX  are  determined  by  using,  in 
addition  to  MO,  two  parameters  DEL  and  IDEL.  DEL,  as  given  by  (4),  is  used  whenever  MO  =  2.  Its 
sign  determines  whether,  for  the  next  cycle,  the  j  index  associated  with  UV  should  be  incremented  or 
whether  the  i  index  associated  with  XY  should  be  incremented.  Also,  when  MO  =  2  IDEL  is  set  to 
DEL/I  DEL  |.  It  is  used  when  MO  =  1  and  its  sign  determines  whether  j  or  i  should  be  incremented  for 
the  next  cycle. 

We  begin  here  with  the  assumption  that  the  first  intersection  point  c  has  been  found  by  FST 
and  the  XY  and  UV  arrays  have  been  reordered  by  SORT  as  described  above.  SOLV  is  called  with 
input  coordinates  yj)  of  vertex  Z|  and  (uj,  Vj)  of  Wj  starting  with  i  =  j  =  1.  Its  outputs  are  MO  =  2, 
the  coordinates  of  the  intersection  point  c,  and  DEL.  We  assume  for  the  first  intersection  point  c  that 
DEL  >  0.  Thus  the  angle  with  vertex  at  c,  measured  in  a  counterclockwise  direction  from  the  line 
segment  cwj  to  the  line  segment  czj,  is  positive.  This  means  that  at  least  a  part  of  ZjSj  wiU  belong  to 
INTX,  so  we  increment  the  index  j.  Thus  in  the  next  cycle  SOLV  will  be  called  to  decide  if  and 
W2W3  have  an  intersection  point.  However  before  the  call  is  made,  INT  is  called  to  see  if  starting  with 
Zj  a  successive  set  of  points  Z  =  {z2,  Z3,  ...,  z„,}  have  already  been  found  by  ITR  (see  below  for  a 
description  of  the  function  INT).  If  this  is  not  the  case  then  SOLV  connders  ZjZ]  and  W2W3  for  the  next 
crossing.  Otherwise  the  set  Z  is  not  empty  for  some  m  such  that  2  <  m  <  NXl  and  i  is  set  to  m.  Then 
SOLV  connders  z„,z„,  ^  i  and  W2W3  for  the  next  intersection.  The  parameter  IDEL  is  set  to  one. 
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If  on  the  other  hand  DEL  <  0,  then  the  roles  of  the  XY  and  UV  edges  are  interchanged  and 
IDEL  =  - 1. 

The  role  of  IDEL  comes  into  play  if  the  output  from  SOLV  at  some  stage  gives  MO  =  1,  which 
implies  that  the  XY  and  UV  edges  under  consideration  do  not  cross.  In  this  case,  if  IDEL  =  1  then  j  is 
incremented  by  one,  and  if  IDEL  =  —  1  then  i  is  incremented  by  one.  No  consideration  is  given  to 
points  previously  found  by  ITR  as  was  done  with  the  Z  sequence  when  MO  =  2. 

If  MO  =  3,  then  SOLV  has  obtained  a  value  for  DEL  that  is  numerically  sero.  Thus  the  two 
edges  under  consideration  are  parallel  and  overlap;  if  they  do  not  overlap  then  MO  =  1  and  XINT 
proceeds  as  described  in  the  preceding  paragraph.  In  case  of  overlap,  the  end  points  of  the  overlap 
belonging  to  INTX  have  already  been  obtained  by  ITR.  The  indices  are  advanced  as  described  for 
MO  =  2. 

FUNCTION  INT-  CALLED  BY  XINT 

The  function  INT(x,  y,  k)  is  used  to  determine  if  x  is  one  of  the  first  k  elements  of  the  WW 
array  and  if  y  is  one  of  the  first  k  elements  of  the  array  ZZ.  If  both  conditions  are  true  then  INT  ^  0; 
otherwise  INT  =  0.  INT  is  called  by  XINT  with  k  =  Kl,  the  number  of  vertices  of  INTX  found  by  ITR, 
(see  page  3),  i.e.,  it  examines  x  and  y  against  the  elements  of  WW  and  ZZ,  which  are  the  arrays 
containing  the  coordinates  of  the  intersection  points  obtained  by  ITR. 
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EXAMPLES  USING  FINTX,  BASED  ON  INTSBC 

Here  we  give  some  examples  of  convex  polygons  XY  and  UV  for  which  the  intersecticm  INTX 
is  determined  using  the  Lahey  Fortran  code,  FINTX,  baaed  on  the  INTSEX7  algorithm.  A  figure  is 
given  for  each  example  showing  the  geometry  of  XY  and  UV.  The  format  of  the  examples  is  explained 
below. 

I,X,Y  specify  the  ith  vertex  and  its  (x,y)  coordinates.  In  Example  2,  the  first  line  refers  to 
x(l)  =  8.0,  y(l)  =  0.0;  J,U,V  specify  the  jth  vertex  and  its  (u,  v)  coordinates.  In  Example  2,  the  sixth 
line  refers  to  u(2)  =  4.0,  v(2)  =  2.0. 

The  next  group  of  data  is  the  result  of  using  ITR.  It  specifies  the  vertices  of  UV  conUuned  in 
XY  followed  by  the  vertices  of  XY  contained  in  UV.  In  Example  2,  the  nineth  line  indicates  that  [u(2), 
v(2)]  or  Wj  is  contained  in  XY  and  is  stored  as  the  first  element  in  WW  and  ZZ.  Note  LW  =  1.  Line  10 
indicates  that  the  first  vertex  of  XY,  Z|,  is  contained  in  UV  and  its  coordinates  are  stored  in  the  second 
element  of  the  arrays  WW  and  ZZ.  Note  LZ  =  1. 

FST  refers  to  the  algorithm  that  determines  the  first  intersection  point  of  INTX  not  found  by 
ITR.  In  Example  2,  P  and  Q  refer  to  the  x  and  y  coordinates  of  this  intersection  point,  which  is  noted 
on  the  figure  by  c,  i.e.,  P  =  x  =  6.50,  Q  =  y  =  2.0. 

Next  MO,I,J,P,Q,DEL  are  given.  They  refer  to  input  to  and  output  from  SOLV.  On  line  14 
of  Example  2,  the  input  is  x(3),  y(3),  u(2),  v(2),  the  output  is  MO  =  2,  the  coordinates  of  the 
intersection  point  x  =  P  =  6.00  and  y  =  Q  =  0.0,  and  DEL  =  32.  Hence  for  Example  2,  SOLV  yields 
the  result  that  the  line  segment  S3Z4  of  XY  and  W3W3  of  UV  have  an  intersection  point  x  =  6.00, 
y  =  0.0.  The  fact  that  DEL  >  0  means  that  J  will  be  increased  by  one  for  the  next  cycle  that  can  be 
seen  in  line  15.  But  I  is  also  incremented  by  one,  since  z^  is  an  element  found  by  ITR  at 
x(4)  =  x(l)  =  8.00,  y(4)  =  y(l)  =  0.0. 

Finally,  INTX  is  given  in  terms  of  the  coordinates  of  its  vertices  that  are  stored  in  the  arrays 
W  and  Z.  For  Example  2,  there  are  four  distinct  points  (vertices  of  INTX),  the  first  of  which  is  stored 
in  W(l)  =  6.00,  Z(l)  =  0.0,  as  shown  on  line  17. 
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EXAMPLE  1 


I,X,Y  1  8.000000  0.000000 
1,X,Y  2  5.000000  4.000000 
I,X,Y  3  0.000000  0.000000 

I, X,Y  4  8.000000  0.000000 

J, U,V  1  8.000000  3.000000 
J,U,V  2  2.000000  .  3.000000 
J,U,V  3  8.000000  1.000000 
J,U,V  4  8.000000  3.000000 

FST,P,Q  5.750000  3.000000 

MO,I,J,P,Q,DEL  2  1  1  5.750000  3.000000-24.0 

MO,I,J.P,Q,DEL  2  2  1  3.750000  3.000000  24.0 

MO,I,J,P,Q,DEL  2  2  2  3.235294  2.588235-34.0 

MO,I,J,P,Q,DEL  1  3  2 

MO,I,J,P,Q,DEL  2  4  2  7.000000  1.333333  18.0 

MO,l,J,P,Q,DEL  1  4  3 

MO,I,J,P,Q,DEL  2  4  4  5.750000  3.000000-24.0 

JJ,W,Z  1  7.000000  1.333333 
JJ,W,Z  2  5.750000  3.000000 
JJ,W,Z  3  3.750000  3.000000 
JJ,W,Z  4  3.235294  2.588235 


EXAMPLE  2 


I,X,Y  1  8.000000  0.000000 
I,X,Y  2  5.000000  4.000000 
I,X,Y  3  0.000000  0.000000 

I, X,Y  4  8.000000  0.000000 

J, U,V  1  11.000000  2.000000 
J,U,  2  4.000000  2.000000 
J,U,V  3  8.000000  -2.000000 
J,U,V  4  11.000000  2.000000 

J,K,WW^Z,LW  2  1  4.000000  2.000000  1 
I,K,WW^,LZ  1  2  8.000000  0.000000  1 

FST,P,Q  6.500000  2.000000 

MO,I,J.P.Q,DEL  2  1  1  6.500000  2.000000-28.0 

MO,I,J,P,Q,DEL  1  2  2 

MO,I,J,P,Q,DEL  2  3  2  6.000000  0.000000  32.0 

MO,I,J,P,Q,DEL  1  4  3 

MO,I,J,P,Q,DEL  2  4  4  6.500000  2.000000-28.0 

JJ,W,Z  1  6UNM000  0.000000 
JJ,W,Z  2  8.000000  0.000000 
JJ,W,Z  3  6.500000  2.000000 
JJ,W,Z  4  4.000000  2.000000 
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EXAMPLES 


I,X,Y  1  10.000000  0.000000 
I,X,Y  2  6.000000  8.000000 
I,X,Y  3  0.000000  0.000000 

I, X,Y  4  10.000000  0.000000 

J. U.V  1  11.000000  5.000000 
J,U,V  2  4.000000  5.000000 
J.U.V  3  0.000000  2.000000 
J,U,V  4  11.000000  5.000000 

J,K,WW,ZZ,LW  2  1  4.000000  5.000000  1 

J,K,WW,ZZ,LW  3  2  9.000000  2.000000  1 

FST,P,Q  7.500000  5.000000 

MO,I,J,P,Q,DEL  2  1  1  7.500000  5.000000-56.0 

MO,I,J,P,Q,DEL  1  2  3 

MO,I,J,P,Q,DEL  1  3  3 

MO,I,J,P,Q,DEL  2  4  3  9.000000  2.000000  28.0 

MO,I,J,P,Q,DEL  2  4  4  7.500000  5.000000-56.0 

JJ,W,Z  1  9.000000  2.000000 
JJ,W,Z  2  7.500000  5.000000 
JJ,W,Z  3  4.000000  5.000000 

EXAMPLE  4 


I,X,Y  1  8.000000  0.000000 
I,X,Y  2  6.000000  6.000000 
I,X,Y  3  0.000000  0.000000 

I, X,Y  4  8.000000  0.000000 

j,u,v  1  ;o.oooooo  4.000000 

J, U,V  2  5.000000  4.000000 
J,U,V  3  0.000000  9.000000 
J,U,V  4  10.000000 

J,K,WW,ZZ,LW  2  1  5.000000  4.000000  1 


J4C,WW,ZZJ.W  3  2  0.000000  0.^30000  0 


FST,P,Q  6.666667  4.000000 


MO,I,l,P,Q,DEL  2 
MO,I,J,P,Q,DEL  2 
MO,I,J«P,Q,DEL  2 
MO,I,J,P,Q,DEL  2 
MO,I,J,P,Q,DEL  2 
JJ,W,Z  1  0.000000 
JJ,W,Z  2  7.058824 
JJ,W,Z  3  6.666667 
JJ,W,Z  4  5.000000 


1  1  6.666667 

2  3  0.000000 

3  3  0.000000 

4  3  7.058824 
4  4  6.666667 
0.000000 
2J23529 
4.000000 
4.000000 


4.000000-30.0 
0.000000  -36.0 
0.000000-32.0 
2J23529  68.0 
4.000000  -30.0 
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EXAMPLES 


I,X,Y  1  8.000000  0.000000 
I,X,Y  2  6.000000  6.000000 
I,X,Y  3  0.000000  0.000000 

I, X,Y  4  8.000000  0.000000 

J. U.V  1  11.000000  4.000000 
J.U.V  2  4.000000  4.000000 
J,U,V  3  2.000000  2.000000 
J,U,V  4  11.000000  4.000000 

J,K,WW,ZZ,LW  2  1  4.000000  4.000000  1 
J,K,mV^Z,LW  3  2  2.000000  2.000000  1 

PST.P.Q  6.666667  4.000000 

MO,l,J,P,Q,DEL  2  1  1  6.666667  4.000000-42.0 

MO,I,J,P,Q,DEL  2  2  3  2.000000  2.000000-42.0 

MO,LJ,P,Q,DEL  1  3  3 

MO,I,J,P,Q,DEL  2  4  3  6.065517  3.103448  58.0 

MO,I,J.P,Q,DEL  2  4  4  6.666667  4.000000-42.0 

JJ,W,Z  1  2.000000  2.000000 
JJ,W,Z  2  6.065517  3.103448 
JJ,W,Z  3  6.666667  4.000000 
JJ.W.Z  4  4.000000  4.000000 


EXAMPLE  6 


I,X,Y  1  10.000000  0.000000 
I,X,Y  2  5.000000  5.000000 
I,X,Y  3  0.000000  0.000000 

I, X,Y  4  10.000000  0.000000 

J. U.V  1  10.000000  2.000000 
J,U,V  2  5.000000  5.000000 
J,U,V  3  2.000000  2.000000 
J.U.Y  4  10.000000  2.000000 

J,K,WW,ZZ,LW  2  1  5.000000  5.000000  0 
J,X,WW^LW  3  2  2.000000  2.000000  1 

FST.P.Q  8.000000  2.000000 

MO,I,J,P,Q,DEL  2  1  1  8.000000  2.000000  40.0 
MO,I,J,P,Q,DEL  2  3  2  5.000000  5.000000  40.0 
M0,I,J,P,Q,DEL  3  3  3  5.000000  5.000000  0.0 
MO,I,J,P,Q,DEL  1  3  4 
M0,I,J,P.Q.DEL  2  4  4  8.000000  24)00000  40.0 

JJ,W,Z  1  2.000000  2.000000 
JJ,W,Z  2  8.000000  2.000000 
JJ,W,Z  3  5.000000  5.000000 
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I.X,Y  1  8.000000  0.000000 
I,X,Y  2  S.000000  5.000000 
I,X,Y  S  0.000000  0.000000 
I.X,Y  4  8.000000  0.000000 


J,U,V  1  0.000000  3.000000 
3,U,V  2  8.OOOO00'  1.000000 
J,U,V  3  3.000000  5.000000 
J.U.V  4  0.000000  3.000000 


J,K,WW^,LW  3  1  3.000000  5.000000  0 
FST.P.Q  6.600007  1433333 


MO,I,J,P.Q,DEL  2  1  1 
M0,LJ.P,Q,DEL  2  2  2 
M0,I,J,P,Q,DEL  2  2  3 
MO,I,J,P,Q,DEL  2  2  4 
M0.I,J.P,Q.DEL  1  3  4 
MO,I,J,P,Q,DEL  2  4  4 


0.000007 

3.000000 

3.000000 

1405217 


1433333  30.0 
5.000000  37.0 
5.000000  9.0 
2.008000  <40.0 


0.000007  1433333  30.0 


JJ,W,Z  1  0.000007  1.333333 
JJ,W,Z  2  3.000000  5.000000 
JJ.W.Z  3  1405217  2.008090 


EXAMPLE  8 


m 


1,X,Y  1  8.000000  0.000000 
I,X,Y  2  4.000000  4.000000 
I,X,Y  3  0.000000  0.000000 
I,X,Y  4  8.000000  0.000000 


J,U,V  1  3.000000  -2.000000 
J,U,V  2  7.000000  3.000000 
3,U,V  3  1.000000  3.000000 
J.U.V  4  3.000000  -2.000000 

FST.P.Q  6.111111  1.888889 


MO,I,J,P,Q,DEL 

MO,I,J,P,Q,DEL 

mo,lj.p.Q.del 

mo,i,j.p.Q,del 

mo,lj,p.q.del 

MO,I,J,P,Q,DEL 

mo,i,j,p,q,i>el 


2  1  1  0.111111 
2  1  2  5.000000 
2  2  2  3.000000 
2  2  3  1471429 
2  3  3  2.200000 
2  3  4  4.000000 
2  4  4  0.111111 


1488889  30.0 
3.000000  -24.0 
3.000000  24.0 
1471429  -28.0 
0.000000  40.0 
0.000000  -40.0 
1488889  30.0 


JJ.W.Z  1  2.200000  0.000000 
JJ,W,Z  2  4.000000  0.000000 
JJ,W,Z  3  0.111111  1488889 
JJ,W,Z  4  5.000000  3.000000 
JJ,W,Z  5  3.000000  3.000000 
JJ.W.Z  0  1471429  1471429 
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EXAMPLES 


I,X,Y  1  0.000000  0.000000 
I,X,Y  2  8.000000  0.000000 
I,X,Y  3  6.000000  4.000000 

I, X,Y  4  0.000000  0.000000  (8.4) 

J. U.V  1  6.000000  -2.000000 
J,U,V  2  6.000000  4.000000 
J,U,V  3  0.000000  0.000000 
J.U.Y  4  6.000000  -2.000000 

J,K,WW.ZZ,LW  2  1  6.000000  4.000000  0 
J,K,WW.ZZ,LW  3  2  0.000000  0.000000  0 

FST.P.Q  6.000000  0.300000 

MO,l,3,P,Q,DEL  2  1  1  6.000000  0.000000-48.0 

MO,I,J,P,Q,DEL  1  2  3 

MO,I,J,P,Q,DEL  2  3  3  0.000000  0.000000-36.0 

MO,I,J,P,Q.DEL  2  4  3  0.000000  0.000000  16.0 

MO,I,J,P,Q,DEL  2  4  4  6.000000  0.000000-48.0 

JJ.W.Z  1  0.000000  0.000000 
JJ,W,Z  2  6.000000  0.000000 
JJ.W.Z  3  6.000000  4.000000 

EXAMPLE  10 


I,X,Y  1  4.000000  1.000000 

LX, Y  2  8.000000  5.000000 
I,X,Y  3  0.000000  5.000000 

I, X,Y  4  4.000000  1.000000 

j,u,v  1  7.000000  0.000000  mja 

J, U,V  2  4.000000  4.000000 
],U,V  3  2.000000  0.000000 
J,U,V  4  7.000000  0.000000 

J,K,WWW>W  2  1  4.000000  4.000000  1 

I,K,WW,ZZ42  1  2  4.000000  1.000000  1 

FST.P.Q  5.285714  2.285714 

MOJ,J,P,Q,DEL  2  1  1  5.285714  2.285714-28.0 

MO, I,J,P,Q,DEL  1  2  2 

MO,I,J,P,Q,DEL  2  3  2  3.000000  2.000000  24.0 

MO.LJ.P.QtDEL  1  4  3 

MO,I,3,P,Q,DEL  2  4  4  5.285714  2.285714-28.0 

JJ,W,Z  1  4.000000  1.000000 
JJ,W,Z  2  5.285714  2.285714 
JJ,W,Z  3  4.000000  4.000000 
JJ,W,Z  4  3.000000  2.000000 
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EXAMPLE  11 


I,X,Y  1  0.000000  0.000000 
I,X,Y  2  5.000000  0.000000 
I.X,Y  3  3.000000  5.000000 
I,X,Y  4  0.000000  3.000000 

I, X,Y  5  0.000000  0.000000 

J. U.V  1  0.000000  -2.000000 
J.U,V  2  5.000000  2.000000 
J,U,V  3  0.000000  0.000000 
J,U,V  4  0.000000  -2.000000 

I,K,WW^,LZ  1  1  0.000000  0.000000  1 
I,K,WW,ZZ,LZ  4  2  0.000000  3.000000  1 

FST.P.Q  2.500000  0.000000 


MO,I,3,P,Q, 

MO,I,J,P,Q, 

MO.LJ.P.Q, 

MO,I,J,P,Q, 

MO,I,J,P,Q, 

MO,I,J,P,Q, 


DEL  2  1 
DEL  2  2 
DEL  2  2 
DEL  2  3 
DEL  2  5 
DEL  2  5 


1  2.500000 

1  4493930 

2  3.823529 

2  2.045455 

3  0.000000 

4  2400000 


0.000000  -20.0 
1.515152  33.0 
2.941176  -17.0 
4.363636  22.0 
0.000000  40.0 
0.000000-20.0 


JJ,W,Z  1  0.000000  0.000000 
JJ,W,Z  2  2.500000  0.000000 
JJ.W.Z  3  4.393939  1415152 
JJ,W,Z  4  3.823529  2.941176 
JJ,W,Z  5  2.045455  4463636 
JJ,W,Z  6  0.000000  3.000000 


EXAMPLE  12 


m 

s 

(0^  ^ 
m  ' 


I,X,Y  1  0.000000  0.000000 
LX,Y  2  5.000000  0.000000 
I,X,Y  3  5.000000  2.000000 
I,X,Y  4  3.000000  5.000000 
I,X,Y  5  0.000000  3.000000 

I, X,Y  6  0.000000  0.000000 

J, U,V  1  0.000000  -2.000000 
J.U.V  2  5.000000  2.000000 
J,U,V  3  0.000000  6.000000 
J,U,V  4  0.000000  -2.000000 

J,K,WW,ZZ,LW  2  1  5.000000  2.000000  0 

I,K,WW^Z,LZ  1  2  0.000000  0.000000  1 
LK,WW,ZZ,LZ  5  3  0.000000  3.000000  1 

FST,P,Q  2.500000  0.000000 


MO,I,J,P,Q,DEL  2  1  1  2.500000  0.000000-20.0 

MO,I,J,P,Q,DEL  2  2  2  5.000000  2.000000-10.0 

MO,I,J,P,Q,DEL  2  3  2  5.000000  2.000000  -7.0 
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MO,I,J,P,Q.DEL  2  4  2  2.045455  4.363036  22.0 

MO,I,J,P,Q,DEL  2  6  3  0.000000  0.000000  40.0 

MO,I,J,P,Q,DEL  2  6  4  2.500000  0.000000-20.0 


JJ,W,Z  1  0.000000  0.000000 
JJ,W,Z  2  2.500000  0.000000 
JJ,W,Z  3  5.000000  2.000000 
JJ,W,Z  4  2.045455  4.363636 
JJ,W,Z  5  0.000000  3.000000 


EXAMPLE  13 


I,X,Y  1  0.000000  0.000000 
I,X,Y  2  6.000000  0.000000 
I,X,Y  3  0.000000  6.000000 

I, X,Y  4  0.000000  0.000000 

J, U,V  1  1.000000  3.000000 
J,U,V  2  2.000000  1.000000 
J.U.Y  3  3.000000  1.000000 
J.U.Y  4  4.000000  2.000000 
J,U,Y  5  3.000000  5.000000 
J,U,Y  6  1.000000  3.000000 

J,K,WW^,LW  1  1  1.000000  3.000000  1 

J,K,WW^Z,LW  2  2  2.000000  1.000000  1 

J,K,WW,ZZ4/W  3  3  3.000000  1.000000  1 

J,K,WW,ZZ,LW  4  4  4.000000  2.000000  1 

FST,P,Q  2.000000  4.000000 

MO,I,J,P,Q,DEL  2  1  1  2.000000  4.000000-24.0 

MO,I,J,P,Q,DEL  1  2  5 

MO,I,J,P,Q,DEL  1  3  5 

MO,I,J,P,Q,DEL  2  4  5  4.000000  2.000000  12.0 

MO,I,J,P,Q,DEL  2  4  6  2.000000  4.000000-24.0 

JJ,W,Z  1  2.000000  1.000000 
JJ,W,Z  2  3.000000  1.000000 
3J,W,Z  3  4.000000  2.000000 
JJ,W,Z  4  2.000000  4.000000 
JJ.W.Z  5  1.000000  3.000000 
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